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ABSTRACT 


A plasma strong turbulence, weak coupling, theory is applied to 
the problem of cosmic ray pitch angle scattering in magnetostatic 
turbulence. The theory used is a rigorous generalization of 
Weinstock's "resonance-broadening" theory and contains no ad hoc 
approximations. A detailed calculation is presented for a model of 
"slab" turbulence with an exponential correlation function. The 
results agree well with numerical simulations. The rigidity 
dependence of the pitch angle scattering coefficient differs from that 
found by previous researchers. The differences result from an 
inadequate treatment of particle trajectories near 90° pitch angle 
in earlier work. 


Subject headings: cosmic rays -hydromagnetics -plasmas 



I . INTRODUCTION 

Central to the problem of understanding the origin of cosmic 
rays is the difficulty of computing the motion of charged particles 
through regions of strongly turbulent electromagnetic fields that 
characterize interstellar and interplanetary space. The basic 
problem is to compute a spatial diffusion coefficient starting 
from knowledge of the statistical properties of the electromagnetic 
turbulence. To do this one needs a valid kinetic theory that 
describes particle motion in strongly turbulent fields. In general, 
such a theory is not available. In this paper, we are interested 
in the special circumstance in which the cosmic ray distribution 
neither modifies the fields with which it interacts, nor undergoes 
Coulomb collisions with other particles. The first partially 
successful efforts to develop such a kinetic theory (Jokipii 1966; 
Roelof 1966; Hall and Sturrock 1967; Hasselmann and Wiberenz 1968) 
made use of the quasilinear approximation (v. , e.g . Rowlands, Shapiro, 
and Shevchenko 1966). In this approximation, it is assumed that the 
deviation of a particle's trajectory from its helical motion in the 
mean magnetic field is small because the strength of the fluctuation 
field is assumed to be weak. It is then possible to compute a pitch 
angle scattering coefficient D^, that describes particle diffusion 
in turbulent magnetic fields (jj. is the cosine of the particles' 
pitch angle with respect to the mean magnetic field). From D one 
can compute a spatial diffusion coefficient (Earl 1974). In this 
paper we will discuss only the derivation of the pitch angle scattering 
coefficient, D^. 
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Recently, several aspects of the quasilinear approximation have 
been questioned (Klimas and Sandri 1973; Jones, Birmingham and Kaiser 
1973c; Kaiser, Jones and Birmingham 1973; Volk 1973). The difficulties with 
quasilinear theory are associated with the behavior of particles with 
pitch angles, 0, nar rr/2, where either D (p = 0) =0 (when the 
fluctuation fields cannot mirror particles to fi.rst order) , or D ((j,=0) 

r 4, 

c x 6 (p.) [when first order mirroring forces do exist (Fisk, Goldstein, 

Klimas, and Sandri 1974; Goldstein, Klimas and Sandri 1975; also 
cf. Lee and VBlk 1975)]. 

In the past two years, several attempts have been made to improve 
quasilinear theory (Jones, Kaiser and Birmingham 1973a, b; VBlk 1973; 

Owens 1974). Owens’ work replaces two of the nonlinear terms that 
are dropped in quasilinear theory with a constant. This results 
in a small amount of scattering at 0 = rr/2, which is insufficient to 
explain the significant scattering through 0 = Tl/2 observed in 
numerical simulations of the problem (Kaiser 1974). VBlk 1973 modified 
quasilinear theory by incorporating a better treatment of the particle 
orbits. He recognized the importance of magnetic mirroring (particle 
trapping) effects in his- description of the scattering. His approach 
is similar to that of Dupree (1966, 1967) for electrostatic turbulence. 

His results are also in rather poor agreement with the numerical simula- 
tions (Kaiser 1974) . This apparently results from a breakdown in 
several of the simplifying approximations made in his analysis.. 
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The most successful theoretical effort to date to improve 
quasilinear theory is that of Jones, Kaiser and Birmingham (1973a, b). 
Their analysis is motivated by the nonlinear (electrostatic) 
turbulence theory of Weinstock (1969). Jones, et al. (1973a, b) 
introduce the concept of partial averaging which enables them to 
include mirroring effects in a physically plausible way. The basic 
approximation made by Jones, et al. ( 1973a, b) is to introduce the 
propagator Up, which describes the motion of particles through 
the partially averaged field. They asrgue, but do not prove, that 

c 

U provides a useful approximation to U . , the exact propagator in 

p A 

Weinstock' s theory. (A mathematical definition of is given in 

§11). The results of these computations appear to agree nicely with 

the three numerical simulations presented by Kaiser (1974). However, 

both the simulations and the theoretical analysis are presently being 

revised (Kaiser, private communication). 

The analysis presented in this paper is a derivation that starts 
with the Vlasov equation and proceeds in a rigorous manner. The 
particle trajectories that are computed are more general than those 
in previous work. We show that these higher order corrections are 
necessary to obtain a good approximation to D at [j. near zero. Our 
analysis follows the ideas of Weinstock ( 1969) , Piran ( 1972) and 
Ben-Israel, Piran, Eviatar and Weinstock (1975). These latter two 
papers are generalizations of Weinstock' s weak coupling electrostatic 
turbulence theory to include electromagnetic effects. As shown 
below, this generalization is not a trivial exercise, and it is not 
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surprizing that the results of VBlk (1973) are quantitatively in- 
accurate. It is not necessary to introduce partial averaging or 
any other heuristic ideas in the analysis that follows. The results 
shown in §V are in good agreement with the simulations that are 
presently available for comparison. 

In the next section we outline the theoretical framework of the 
turbulence theory under the fundamental assumption that the coupling 
is weak (Kadomstev 1965). Much of the general theory has been developed 

by Piran (1972) and Ben-Israel, ej; al. (1975) and we will not repeat 

* 

their detailed formulation here. The general theory leads to equations 
that cannot be solved without limiting one's attention to specific 
models of wave turbulence. One such model of interest for the problem 
of cosmic ray pitch angle scattering, is homogeneous "slab" magnetostatic 
turbulence with an exponential correlation function. s 

Because cosmic ray velocities are typically much greater than the 
phase velocities of the magnetohydrodynamic waves that scatter them, 
the magnetostatic approximation should be very well satisfied. This 
is equivalent to dropping the Maxwell "curl" equations in the analysis. 

In confining our attention to "slab" turbulence, we consider only 
magnetic fluctuations that vary along the direction of the mean magnetic 
field. This greatly simplifies the analysis. Waves observed in the 
interplanetary medium near 1 AU are thought to be predominantly Alfv&n 
waves propagating along the mean field direction. A magnetostatic slab 
model approximates this situation quite well. The specific examples 
computed below are limited to an exponential correlation function for 
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the turbulence. This is done to compare this theory with numerical 
simulations of this problem that have thus’ far been confined to an 
exponential correlation function (Kaiser 1974)., Generalization to 
other, correlation functions and to isotropic turbulence is straight- 
forward and will be considered in future papers. 

Section III is devoted to a detailed calculation of D (|j = o) 
valid to second order in the strength of the fluctuation field* All 
nonlinear terms are included except those that are proportional to the 
fluctuation in the parallel mementum. They are expected to be un- 
important for the reasons discussed below. 

In section IV, we estimate D at all values of m. In doing 
so, we neglect most terms that are proportional to the deviations of the 
mean of the momentum from the value expected from quasil inear theory. 
The validity of this approximation is estimated by comparing the 
results of section IV at (j. = o with those of section III. For the 
range of parameters considered, the approximation is good. 

In the last section, we compare our results to the simulations 
and to previous theoretical work. We find that the results of 
VBlk (1973) and Owens (1974) fit neither the numerical simulations 
(Kaiser 1974) nor our theory. The published numerical simulations 
do not cover a sufficiently wide range of parameters to distinguish 
between this theory and that of Jones, et al. (1973a, b). However, 
the two theories predict different dependencies on the parameters 
of the theory. We show that the results of Jones et _al. (1973a, b) 
can be recovered from our analysis if we drop terms important for 
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particle mirroring. We conclude that our analysis contains the first 
accurate evaluation of particle orbits near 90° pitch angle. 

II THEORY 

The derivation presented here follows, in broad outline, the 
work of Piran (1972) and Ben-Israel, Piran, Eviatar and Weinstock 
(1975). We begin with the Vlasov equation, written in dimensionless 
form (Klimas and Sandri 1971, 1973), 


— + KF + £F + T)£ / J? = 0 ( 1) 

dT 

for the cosmic ray distribution function F(x,£,t). In equation (1) 

T *» to) where t is dimensional time and tw is the gyrofrequency of 
c c 

the particle. In addition, the differential operators K, £ and £* 
are defined as follows: 


K 


. p.3 


dx 


a a 




( 2 ) 

(3) 

(4) 


where £ = £/p,x = x /r (x* is dimensional length and r is the 
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Larmor radius) and 



The skew -symmetric tensors Q and Q* are defined as 0^ 

and O' - € S' where £ = B /B , 0* = B' /<B' *B'> 1/2 and the 
ij ijk k “o o 

magnetic field is separated into a mean part, j^, plus a fluctuation, 
B'. The small parameter is T| = -B'> 1/,2 /B . The symbol <> denotes 

— . — o 

ensemble average. 

We proceed to derive a diffusion-like equation for the average 
part of the distribution function f = <E>. We assume that Kf ~ 0. If 
one takes the ensemble average of equation (1) and subtracts the 
resulting equation from equation (1), the following relation for 
J: ' = F - results: 

[§f + (1 - f ' = "^ £,f (6) 

where W = K + £ + T)£' is the total Hamiltonian differential operator 

and the operator A takes the ensemble average of everything to its 

right (i.e. AF = O) . Owens (1974) takes the nonlinear terms (1-A) 

TjX’ on the left hand side of equation (6) and sets thon equal to a 

constant, a. This severe approximation results in a gross undirestima 

tion of the scattering through (j. = o. 

Equation (6) can be solved formally in terms of the propagator 

U (T, t ), which is defined by the following differential equation 
A o 

(Weinstock 1970; Ben-Israel, et al. 1975): 

+ (1 - A)>0 U A (T - T„) = 0 (7) 

U (0) = 1 
A 
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In writing equation (7), we used the fact that for magnetostatic 

turbulence U A (t»T ) = U A (T-t ). 

** o A o 

Equation (6) has the formal solution 

T 

f'(T) = -T] J d\ U A (t-O £’f (\) (8) 

T Q 

where we have taken f'(o)=0. 

If we substitute this solution back into the ensemble average of 
equation (1) , we have the exact master equation for f: 

+ if = T] 2 f dX U A (T-X) £<> f(X) (9) 

T o 

Because it is not possible to find an exact representation for 
U A <T-T 0 ), one proceeds by finding approximations that retain as much 
information as possible about the effects of the stochastic fields 
on the particle trajectories. With this goal in mind, we do not 
make the quasilinear approximation in which it is assumed that the 
stochastic fields do not affect a particle's motion during an inter- 
action. 

Note that equation (9) can be rewritten as 

T 

+ if = if < r r d\ u. (t-x) (i-a) r> f<\> do) 

dT A 

T o 

An approximation to (t) (1“A) can be found in terms of the 
propagator U(t-t q ) defined by 
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( 11 ) 


[&- +X] U(t-T.) = 0 
St 0 

U(0) = 1 * 

This can be written (Ben-Israel, et al. 1975): 

t§- + (1-AjV] (l-A) U(t) = -T|( 1-A) £'AU(T) 


with the formal solution 


T 

(1-A) U(t) = U a (t) (1-A) -“H J d\ U a (t-\)(1-A) X'AU(X) 

o 


( 12 ) 


(13) 


Here, and in the remainder of this paper, we take T q = 0. To proceed, 
we drop the second term on the right-hand side of equation (13). This ’ 
approximation is valid in the limit of weak coupling (Kadomtsev 1965; 

Klimas and Sandri 1973), i. e. in the limit T) « 1. The representation 
of U( T ) is also unknown, so it is necessary to approximate U(t) in terms of the 
propagator U = <U > , which satisfies the ensemble average of equation 
( 11 ): 

If D(T) = - (14) 

or 

C|- +X 3 U(t) = -t] < £'U(t)> + 1} (15) 

OT 

The formal solution is 

T 

U(T) = U(T) + T1 / dx U(t-X) [ X’U(X) - <£'U(X)> ] (16) 

o 
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Therefore, to lowest order in T],U(t) U(t) and equation (10) can 
be rewritten as 

+ £f ~ n 2 < £' / dx U(\) £’> f(r-X) (17) 

o 


To reduce equation (17) to a diffusion equation, we make the adiabatic 
approximation. The nature of this approximation is to restrict 
our interest in equation (17) to times T » 1 when f(T-\) =* f(T). 

One now substitutes equations (3) -(5) into equation (17) and 
averages over gyrophase (c£. Goldstein, Klimas and Sandri 1973). 

After some manipulation, equation (17) becomes 


Mill 

St 


S_ 

an 


D f- f(T) 


(18) 


where 


D (T) 

Id 




^0 <p-Q*i'(x) u(x) g.'(x)-n*p> 


(19) 


D = Lim D (t) 

Id Id 


( 20 ) 


T “* oo 

In writing equation (18), we assumed that the turbulence is axially symmetric 
about f§. Also f ( t) in equation (18) now denotes the gyrophase average 
of f ( t) in equation (17) . 

The assumption in writing equation (18) is that equation (19) 

converges 'rapidly" to D . However, in quasilinear theory, for 

M* 

particle pitch angles near jt/2 , equation (19) takes an arbitrarily 
long time to saturate. (For a more, comprehensive discussion of the 
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adiabatic approximation, we refer the reader. to Weinstock 1970; 

Jokipii 1971; Klimas and Sandri 1973; Jones and Birmingham 1974; and 
Fisk, at al. 1974) . A consequence of the present nonlinear analysis 
is that the integral in equation (19) saturates "quickly 11 , even for 
(j = 0. Thus we expect the adiabatic approximation to be a good one. 
This will be examined in more detail in a future paper. The remainder 
of this paper will be devoted to the derivation of D^(t)» 

We proceed by writing j3'(x) as a Fourier integral transform over 
wave number, k. It folldtos that 

<i'(k) i'(k')u exp[ik'*x ]*Q*p > = 

( 21 ) 

< i.’(k) g.'(k') > <exp [ik'-x(X)] *Q*p (\) > 

This permits writing the theory solely in terms 
of the two-point correlation tensor of the magnetic turbulence. 


We define R(k) through the relation 

3/2 

I'(k') > = (2tt) R(k) 6 (k+k 1 ) (22) 

where R(k) is the Fourier integral transform of the two-point correlation 
tensor, 

R(.l) = < i’<x + r) > (23) 
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(In writing equation (22), we have restricted ourselves to homogeneous 

turbulence.) D (t) can now be written as 

M 1 


-n 2 T 2n 

d u (t) = 2 r r r exp (*£*£) p* 

^ (2tt) “ : 


P(\) exp C-ik*x(\)] > 


(24) 


The basic mathematical problem is that of computing the ensemble 
average in equation (24). If £(t) is written as the sum of a fluctuation 
part, p'(t), plus a mean part < p(t) >» then the correlation becomes 


< p(t) exp [-ik*x(T)] > = < p(t) > < exp C-ik*x(T)] > 


+ < p(t) exp[-ik*x'(T)]> 


(25) 


We define 

A£(\) = < p(\) > ~ P Q (X) (26) 

where p (\) is the streamed velocity of the particle in the absence 
o 

of fluctuation fields. The propagator of this trajectory, IT(t)» 
satisfies 

[-^ + £ ] U 0 <T-T 0 ) = 0 (27) 

U (0) =1 
o 

The solution, is the well-known helical trajectory used in quasilinear 
theory, which we can write as (v, , e.g . , Goldstein, Klimas and Sandri 
1975) . 
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( 28 ) 


U q (t) x a x o (t) = x-r Q (T) 
r Q (T> (t)*P a p* T) (29) 

U q (t) P = P q (t) = <J (T)*p = p*($ (t) (30) 

C (t) = P + N cos T * 0 sin T (31) 

<0 (t) » Pt + N sin T + Q (cos t-1) (32) 



(33) 


N. . 
3-J 



(34) 


Equation (27) can be rewritten (Ben-Israel, et al. 1975) as 

[ l^+V] U 0 (T) = 1) £*D 0 (T) (35) 

with the solution 

T 

U o (t) = U(t) + 11 J* d\ U(T - A) X'U 0 (A) (36) 

o 

The quantity A£>(\) in equation (26) is unique to the magnetic 
turbulence theory. In the electrostatic problem (Dupree 1966; 
Weinstock 1969, 1970), the momentum trajectories do not appear 
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explicitly. The velocity dependent nature of the magnetic force 
greatly complicates the computation of the nonlinear wave-particle 
interaction in strong turbulence (Ben-Israel, et al. 1975). However, 
some simplif ications are possible for the test particle problem in 
magnetostatic turbulence. For example, we will argue below that 
the second term on the right-hand side of equation (25) is unimportant 
in the particular examples that one can compute. We now outline the 
computation of < exp -i [ _k* xCX) 3 ^ and Aj>(X) . 

Following Ben-Israel, et al. (1975) and Weinstock (1970) , we 
expand ■< exp -i [k*x(X) ] > in a cumulant expansion (Kubo 1962). To 
second order we have 

< exp-i [k*x(X)3 > = exp [-ik*< x(\) > - ~ k’Dj^x) *k] (37) 

where D X (\A) = D^X) s < x'(\) x'(X) > and x(X) = (1"A) x(?0 • [Ben-Israel, 
et al, (1975) are missing the factor of 1/2 on the rhs of equation (37)], 


To our knowledge, the general convergence properties of this 
expansion are unknown. However, if the fluctuation, x'(X), obey 
Gaussian statistics, then equation (37) is exact. 

We define 


Ax(T) = < x(t) > - X (t) 

where, again, subscript M o” denotes streaming by the propagator 
U q (t). Then, from equations (36), (3) (4) and (16), we find to 

<7oi 2 > 


(38) 


Ax(t) = -T| 2 J 1 dx J* dv < U(t-v)p*0' 

° \ + 
|tU(v-x) p-n'-g (\) >■ 


(39) 


* 
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Similarly, 


T T p, + 

A£(t) = T] 2 J* d\ J dv < U(T"v)p-Q'* ^ u<v-\)p •$' *Q (X) > 


(40) 


To confute D^(t) , note that one can write 

U(t) = U q (t) - T] fd\ U o (T-\) < £'U(\) > - U q (T) + O (ll 2 ) (41) 

o 

Then, from the definition of D^(t) and equation (41), one has 


D^T) - Tl 2 < J* dv 


where the symbol 


u(t-v) p*n' '£f (v) 


J d\ u(t-x) p -d-£f (x) > 


i 


(42) 


means that operators within the bracket act 

only on variables contained within the bracket. 

To proceed, and at the same time keep the mathematics tractable, 

we will make several of the simplifying assumptions alluded to in the 

introduction. The interested reader is referred to Ben-Israel, et al . 

(1975) for a more general derivation. However, Ben-Israel et al. 

* 

make an unnecessary approximation, viz . , they write 


2- p(T) = [ < p(T) > + p'(T)] f-r P C (T) (43) 

The term “T Ap(r) > neglected in writing equation (43), is of the same 
op - 

order as A£(t), which is kept in their analysis. We will find below 
that an approximate solution can be found, including terms proportional 
to dAp(T)/dp f for the special case of Q = tt/ 2 (|j. = 0) in a "slab" model. 
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We will follow Ben-Israel et al. (1975) in dropping the term 
proportional to ^7 p(t)» We confine the remainder of the discussion 
to the "slab" correlation function 


R(k) = N 6(k l )R(k || ) (44) 

where R(k ) is the power spectrum of the turbulence and 


Jdk^ 6(kj_) = 2 tt 


(45) 


In the following section, we develop the theory for the single 
point |j, = 0. This is the region in which quasilinear theory gives 
the poorest approximation. It is also the case for which the most 
complete nonlinear analysis can be done. 
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III. SCATTERING THROUGH y,*0 IN A SLAB MODEL 

An immediate consequence of restricting the discussion to p,=0 in 

a slab model is that Ax (t) does not appear, Ax (t) = 0 and the second 

l t! 

term in equation (25) can be shown to be proportional to (j, and, hence, 
equal to zero. Using equations (24)-(26), (37) and (44)-(45), Dpt ( t) 
can be written as 


p . r 2 tt » 2 

D(t) = T1 JdX Pd0 Jdk R(k )[cosX + p *Ap (X)]exp[-k D.. (X.) /2] 

* 3/2 o o -co .11 'I l-l II x 


(2tt )' 


(46) 


where 

D 1 (X) = D 1 (X). 

The two quantities that must be computed are Ap (t) and D 1 (t) . 

“1 1 

We outline first the derivation of D^(t) to illustrate the techniques 
used in evaluating all of the nonlinear quantities presented in this 
discussion. We start with equation (42). In general, Q' can be 
expressed as 

O' = ($•£')£ - £(£'.£) - (9-0P. (47) 


In the case of slab turbulence considered here, the first term in 
equation (47) is zero. D-^(t) becomes 
T T 

D l(T) = 'll 2 jdv jdX vX Cu(t-v) P-Q-&* H U(t-v) p-0.£'> (48) 

o o 

One now writes the fluctuation fields as Fourier integral transforms, 
and uses equations (21) and (22) to arrive at 


d x (t) 


T T 

jdv jdX Jd 3 k -<exp[i k* (x(t- v)-x(t-X)) ] vX 
(2tt) 3/2 0 ° 

p(T-v) *fl*R(k) .fj.p (t-X)> 


(49) 
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Substitution of R(k) from equation (44) yields 
-2 f T 00 

d i ( t) = i* dv i dx J dk n R(k it ) vx ^expClkuxJ (v-X)]p^(T-v) •p 1 (T-X)> 

O O "OD 

(50) 

where we have used the approximation (Ben-Israel, et al. 1975) 

Xj'Ct-v) - x(|(T“X) ~ x[|(v-X) (51) 

To evaluate the correlation that appears here, one must expand the 
integrand of equation (50) as follows 

! '«<exp[ik||X 1 , i (v-X)]p^(T-v)*p ± (T-X)> = ■<p^(f-v)> i *<P 1 (T-X> x 

^expik^x* ( V -X)> -f- -<p^(t-v) *P^(t-X)> ^expCik^xf, (v-X)] + 

<pl< T -v) e x p [i k ! x' ( v -X)]>^<p^( T -\)> + 

<p^(T-X)exp[ik j! x j ' 1 (v-X)]>i<p i (T-v)> (52) 

In Appendix A, we show that for |j,=0, only the first term in equation 
(52) is expected to be important. Consequently, equation (50) becomes 
2 T T 00 

D 1 (t) = -2— jdv jdX Jdk R(k )vX expC-k^ (v-X) ][p (t-v)+Ap (t-v) ] 

/5tt o o -od 11 " 1 1 -1 

*[p 01 (T“X) + Ap ± (T-X)] (53) 

Finally, we must derive an expression for Ap ( T ). From the definition 

p 

of Ap^(t) (equation 26) and the identies (35) and (16), we have to ©'(Tj ) 

T T 

ip .00 = T| 2 rax jdv ; -U(T-v)p-CJ'-S- U(v-X)p-n' p l0 (X)v (54) 

o X dp 3p 

In Appendix B, we show that equation (54) may be written approximately as 

- l(T) = % I dX 1: ip,i< T - x )ll= 0 

[p ij0 ( T " v ) + ^P j^T-v)] ’C.* (\) (55) 
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The approximations involved in deriving equation (55) are discussed in 
Appendix B. In Appendix C, a derivation of Ap^r) is given. From that 
discussion we have 

2 t 00 

- Ap (t)| " ^ JdX Jdk j! R(k il )cos\(T-X)exp[-k^D 1 (X)/2] 

** |i,=0 y^TT O 


*{2+kfM> - }| 


(56) 


ll,=0 


where BAx^t)/^ is given by equation (D28) . 

The set of equations (53), (55) and (56) provides a complete system 
of coupled, nonlinear integral equations, which in principle can be 

solved, and substituted into equation (46) for D (t). Unfortunately, 

M- 

this system cannot be solved analytically. To solve the system 
numerically, we must specialize to particular forms of the power spectrum, 
R(kjj). The numerical simulations developed by Kaiser (1974) are for 
R(kjj) of the form 


R (k it ) = (1 + ^k 2 )“ 1 

" JTv 11 


(57) 


where € is the ratio of the correlation length of the magnetic field 
turbulence, \ , to r . If we use equation (57), the integrals over k.. 

C § 11 

can be performed analytically. For these two reasons (comparison with 

numerical simulations and mathematical simplicity), we have confined our 

numerical analysis to the particular case of equation (57) . There are two 

integrations over k- to perform. They are 

2, 

j, u;/ 2 

(58) 

• d i(t) 1/2 


g (t) = J’dk lj R(k 11 )exp-k^D 1 (T)/2 
■00 

= /2t\ Erfc exp[D 1 (T)/2€ 2 ] 
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and 


h ( t ) = Jdkjk 2 * R(k Jexpr-k^D^T)/^] 

0 (59) 

= (2tt€ 4 )” 1/2 (1[2 Tt€2 /d 1 (t)] 1/2 - TtErfc [ D l< T >/a^] 1/2 exp[D 1 <T>/aS 2 ]} 
The system of equations (46), (53), (55) and (56) can now be simplified 

-2 T 

D (t) = J — q/p to J*dXg(X) [cos\ + cos0 Ap (T)+sin0 Ap (t)] (60) 

P (2n) J/£i o o x y 


«f|2 T T - . 

D- (t) = -* — to JdX g(v-X)vX[p 0 (T-v) + Ap (t-v)]*[p 0 (t-X)+4>(t-X)] 
1 O O 1 -*■ 1° -x 

(61) 

2 T T 

Ap (t) = ^ JdX jd v g(v-X) {l+[-- Ap (t-X)] -i x 
- 1 ySo i ^ 11 ^L=O j 


[^ i0 (t-v) + Ap^(t-v) ] -C* (X) 


(62) 


.2 T 


^7 Ap Or) | ='^ IL jdX(T-X) cosX {g(X) + Xh(X)[X- |-Ax (X) ]} (63) 

|i=0 y^TT O ^ |Jb=0 


2 T 

Ax ( T )| = JdX cosX(T-X) 2 {g(X) + Xh(X) [X - ■“ Ax (X)]} (64) 

11 (j,=o rf o " p,=o 

Several additional simplifications can be made. For the axisymmetric 
turbulence considered here, the integrand of D (t) should be independent 
of 0. From the form of equations (60)-(64), this assertion may not be 
obvious. However, it can be shown numerically that this is indeed the 
case-. Therefore, to simplify the analysis presented here, we use the 
fact that the integrand of equation (60) is independent of 0 and 
evaluate it at 0 =0. The following time symmetries are then not difficult 
to prove * 
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d^t) = d 1 (-t) 


Ap (t) = Ap (-t) 

X A 


Ap (t) = -Ap (-t) 

y y 


i Ax ii (T) ■ - * Ax n ( - T) 


(65) 

( 66 ) 

(67) 

( 68 ) 
(69) 


Equations (60) -(62) simplify to 

T 

fax g(X) [cosX + bp (X)] 
JL X 


d.(t) = X— 

/2tt o 


M- 


(70) 


2 <r 

D,(t) = JdX cosX g(\) [t 3 /3 - X/2(t 2 -X 2 /3)] + 
/2rr o 


2 T T 

j d ^ ,f dX S(v-X) vX {Ap (t-v)[2cos(t-X) + Ap (t-X)] + 




o o 


± x 


Ap^y(T-v)^sin(T-X)+ Ap^ y (T-X)]} (71) 


/AP iX (f)v 2 T T 

( ) = fdX fdv g(v-X) [1+ — Ap (t-X)] 

\a P (t)/ fl TT j O 1 ^ 11 ^ C 

±y 


!(. 


cos(t-v+X)\ / Ap v (t-v) > 

1+cosX I 

sin(T-v+X)/ \Ap 

±y 


> iX (T-v)\ 

>„(T-V)/ 


+sinXl 


/ A P iy < T -^\) 

'-Ap v (t-v)// 


(72) 
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To produce a numerical solution, one would like to transform 
these equations into a system of differential equations. However, the 
best that can be done with this system is to derive a Set of nonlinear 
integro- differential equations which can, however, be numerically 
integrated using a standard Adams-Bashforth predictor-corrector technique 
(Krogh 1966). In Appendix D, we outline the transformation of the system 
of equations (70)-(72) and (63) - (64) to a set of integro- differential 
equations. We defer until §V a discussion of the numerical results of 
that integration. We turn now to an appropriate derivation of the system 
of equations that describes scattering at arbitrary pitch angle. 
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IV SCATTERING THROUGH ARBITRARY PITCH ANGLE - AN APPROXIMATE 
CALCULATION IN A SLAB MODEL 


In this section, we derive a set of approximate equations to 
describe pitch angle scattering at arbitrary pitch angles in a slab 
model. Many terms proportional to Ap (t) will be dropped in this analysis 
for reasons that are basically heuristic. An estimate of the validity of 
the various approximations made below can be made _a posteriori by com- 
paring the results of §111 with the results of this section. The system 
of equations derived here is similar to that derived by Ben- Israel, et al . 
(1975) in their discussion of ion-cyclotron turbulence. 

The equation for D (t) is now 


V T) 




2tt 00 


(2TT) 3 / 2 o o 


JdX Jd0 Jdk R(k ) cos [k (p,X-Ax (X) ] exp[-k 2 D X (X)/2] 


[(1-|jl 2 )cos\ + Ap (X) • p ] 

i. i 


(73) 


The first quantity we evaluate is Ax^t), defined by equations 
(38) and (39) . When one uses the approximation (43) , Ax,, (t) becomes 


Ax„(t) = 


jdX J*dv Jd~*k JcT^k' <expi[k»x(T-v) *f k- x(t-X) ] 


(2tt) o x 


(74) 


P Q (t-v) .Q’ (k) • [£ + (v-X) -O' (k' ) + i/j^v-X) -kp 0 (t-v) £ + (v-X) (k' ) ] -$X> 

In writing equation (74), we have replaced p(t~v) by p Q (T-v) 
everywhere. The justification for doing this is two- fold. First, one 
expects that nonlinear terms such as A£(t) and p'(t) will be less 
important for p, not near zero. The reasoning here is that quasilinear 
theory should be a very good approximation except near p, = 0. As we 
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shall see in §V, the results of this section are generally in good 

agreement at p, = 0 with those of §111. Second, the term Ax^(t) is 

equivalent to a shift' in the wave-particle resonance (Weinstock 1972; 

Ben- Israel, et al. 1975), which is generally not the most important 

term in resonant broadening theories. The most important term 

characteristically is the broadening of the wave-particle resonance, 

which saturates the interaction and allows particles to propagate 

through p, = 0. The resonance broadening is given by the term 
2 

exp- [k || D l (t)]. The contributions to D^(t) from Ap(t) and p' ( t) 
should be even less important for p, ^ 0, where quasil inear theory becomes 
a gt/od approximation. 

One now substitutes the identity (47) into equation (74) and uses 
equations (22) and (44) for slab turbulence. After some rather lengthy 
matrix algebra, Axjj (t) becomes 
2 T T oo 

Ax |; (t) = JdX jdv Jdk,| exp{i k^[n(v-X) - Ax )| (v-k)]-k (| 2 D 1 (v-X)/2} x 
^(kjpXcos (v-X) [2p,-i k n (v-X)(l-n 2 )]} (75) 


Because the integrand in equation (75) is basically a function of 
(v-X), AX||(t) can be further simplified to 


2 T 


Ax || (t) - JdX Jdk R(k ) exp[-k ^D 1 (X) ]cosX( T -X)' 
/2 tt o 


(76) 


{2pjCos£k|| (p,X - Ax || (X)) ] + (l-p i 2 )k| ) Xsin[k| 1 (^X - Ax^X)) ].} 


Piran (1972) and Ben-Israel, et al. (1975) incorrectly write the 
argument of the trigonometric functions as [pA+Ax^X)]. 
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W&en R(k||) is given by equation (57), the integrals over can 
be performed analytically with the result 
2 T 

Ax. Or) = jdX(T-X) 2 [2|iG(X) + (1-u, 2 ) H(\)] (77) 

11 /2tt o 


where 

00 

G(t) = Jdk R(k)cosT cos[k((j,T-Ax (T>)]exp[-k 2 D 1 (T>/2] (78) 

n '' 


= cos-r e 2 {e y Erfc[/z -y/a/z]+e y Erfc[/z + y/a/z]} (79) 

a/^TT 


and 

H(t) = JdkR(k)krcosT sin[k(p/r - Ax (T))]exp[-k 2 D 1 (T)/2] (80) 

o 


= 2 ^/^ ^ cosT exp(z^){e" y Erfc[z-y/2z] -e y Erfc[z+y/2z]} (81) 

with 

z s D x (t)/^€ 2 (82) 

y = [|J»T - Ax )( (t)]/€ (83) 


To evaluate (t) , one substitutes equation (44) for R(k) into 
equation (49) . After some algebra, one finds 


,2 T 


D 1 (t)= Jdv J'dk.jvX(l-|i 2 )R(k )cos (v-X) exp[i k [|i(v-X) - 
x /2tt o o -oo 1 


Ax (1 (v-M] - k * l D 1 (v^X)/2] (84) 

where we have again set p(r) = P q (t). 

The double integral over time can be reduced to a single integral 
so that D^(t) becomes 
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D i(t) = JOL (l- |i) 2 ){| T 3 JdXG(X) + JdXG(X)[ T 2 X - x3/ 3] 

/5tt -t 


-T 


- jdXG(X)[T 2 X - X 3 /3]} (85) 

o 

By differentiating D^(t) with respect to t, and setting t' - -T, one 
can easily show that D,(t) = D^(t' - -t) • Consequently, 

D. <t) « ^ (l-|j, 2 ) fdXG(X) [§ t 3 -t 2 X + x3 /3] (86) 

1 JTu J o 3 

The reader may have noticed that although we have dropped Ap (r) 

in our derivations of Ax^(t) and D^(t), we have kept a term 

proportional to Ap (t) *P in the expression for D (t) (equation 73), 

X i M 1 

We expect that deviations of a particles momentum from the predictions 

of quasilinear theory would be relatively unimportant in computing 

the resonance shift (Ax^(t)) or the resonance broadening exp-[k^ D^(t)}« 

However, in computing D (t) a contribution from Ap(t) should be important 

M* — 

for the following reasons: in quasilinear theory for p, = 0, the 

wave - particle resonance is with waves of zero wavelength, which have 

zero power. Therefore, D = 0 at (j, = 0 for the slab model in quasi- 

P 1 

linear theory. The functions D^t) and Ax (t) serve to broaden this 

resonance and ship, it away from k = «. This leads to a finite value 

for D at ^ = 0. In addition to the resonant effects (k^ ~ co) particles 

P 1 

are affected by long wavelength fluctuations (k^ ^ 0) that tend to mirror 
them (Volk 1973; Fisk, et ail. 1974; Goldstein, Klimas and Sandri 1975; 

O 

Lee and Volk 1975). For k^ ~ 0, neither Ax (t) nor exp[-k^ D^(t) ] 
provides significant contributions to D . The only large nonlinear 


27 



contribution to for k^ ~ 0 is the term proportional to Ap^(t) in 
equation (73). This term was dropped by Ben-Israel, et al. (1975) and 
similar terms were neglected by V81k (1973) and by Jones, et al. 
(1973a, b). We will return to this point in §V. We can approximate 
this term in much the same way in which Ax^(t) was computed above. If 
one drops all factors of p'(t) and Ap(t) in the integrand of equation 
(55), one can easily show that 


A p 

Ap Or) *p = (1-m- )cost6p (t) 
1 ^ 


(87) 


2jf' p T 

(l-|i ) cost JdX(T-X)cosX(l+tanX tanT) [G(X)-|j,H(X)‘] 

v/2tt o 


so that 


2 2 T 

D(t) ° ) jd\G(X)[Wp (\)] (88) 

^ /StT O 

The equations in the system (77)-(83), (86)-(88) now form a 

coupled nonlinear system of integral equations for D (t) . Volk (1973) 

derives an expression similar to equation (88), but he neglects the 

term equivalent to 6p (X) and sets D (M<) equal to a constant when 

1 P- 

evaluating «<exp-i k l| x |) (T)> (equation 37). He then finds Ax (t)= 0 and 

1 3 3 

D^(t) - ^ ® T . We find that D, (t) t only when t » 1, and consequently 
the magnitude of D 1 (t) is generally much smaller than would be predicted 
by Volk (1973). In general, his results do not agree with ours, nor do 
they fit the simulations (Kaiser 1974). 

In Appendix E this system of integro-differential equations is 
rewritten as a system of coupled differential equations that can be 
integrated numerically using standard techniques. The results of that 
integration are discussed in the next section. 
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V. DISCUSSION 


The equations that determine D were integrated numerically and 
the results are shown in figures 1-5. For comparison, we have included 
in figures 1-3 the theoretical prediction of Jones, Kaiser and Birmingham. 
(1973a, b), V51k (1973) and Owens (1974). In all of these figures, we 

have plotted £D to conform with the normalization used in the simulations. 

M* 

[If one were to define D to be the dimensioned value of D , then 

£D s (X /uj r )D .1 The theoretical predictions of quasilinear theory 
|ju c c g jj, J 

(Jokipii 1971) are also shown in figures 1-3. The simulation results 
are taken from Kaiser (1974). Since that reference was written, the 
simulation program has been subject to some revision (Kaiser, private 
communication) and the simulation results plotted in figures 1-3 may be 
slightly inaccurate. Consequently, in our discussion, deviations of 
0% of our theoretical curves for €D vs. (j, from the simulation points 
should not be considered significant. 

In figures 1-3, the solid squares are the results of integrating 
the complete set of equations derived in §111 and Appendix D. The 
solid curve results from integrating the system of equations derived in 
§IV and Appendix E. 

The factor of /2 in front of T] in the figures arises from the 
difference between linearly polarization turbulence (Jones, et al . 

1973a, b) and the planar polarization used here. The parameter of 
Jones, et al. (1973a, b) is equal to /2T]. However this equivalence 
between our axisymmetric model and the linearly polarized model used 
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by Kaiser (1974) in the numerical simulations cannot be exact. 

Differences arise when the phase average is performed on equation (17). 

In general, for nonaxisymmetric turbulence, there would be an additional 
term in equation (18) of the form D^-^- f(T). For axisymmetric 
turbulence is not a function of 0 and equation (18) follows. [One 
cannot avoid this complication by assuming that f is gyrotropioij Jones 
(private communication) maintains that the consequences of this lack 
of axisymmetry in the numerical simulations is expected to be small. 

We proceed below to compare our results with the simulations on the 
assumption that this is indeed the case. 

The results shown in Figures 1-3 indicate that quasilinear 
theory is quite good for \x > 1/2, as expected. Significant 
deviations from quasilinear theory are evident for jj, ~ 0. Here our 
results are similar to those of Jones et al, (1973a, b) and Kaiser (1974). 
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The predictions of Volk (1973) are about a factor of two too high for 
/2T| •= .3 and 6 = 1. The' results of Owens (1974) are consistently low 
and will not be discussed any further here. All of these theories have 
very different functional dependencies of the parameters T| and € at jj, = 0. 

Figures 4 and 5 show plots of €D (p, = 0) vs. 6 and 1]. One can see that 

M* 

€D oc T] 3 € (89) 

M- 

whereas, from the summary of Kaiser (1974) one has 

0> JKB ccT1 3 (90) 

M- 

€D V5LK cc Tl 4 /€ (91) 

P* 

€D 0WENS oc €7} 4 /(l+T| 4 ) (92) 

M- 

The parameter dependence illustrated in equations (91) and (92) 
results from the severe approximations made by V*dlk (1973) and Owens 
(1974). Within the context of the turbulence theory presented here, 
there appears to be no consistent way to recover the results of Volk 
(1973) and Owens (1974). However, a detailed comparison between this 
work and that of Jones, et al. (1973a, b) is possible. One can recover 
the parameter dependence shown in equation (90) from our theory when 
one uses in our analysis the approximations to the particle orbits used 
by Jones, at al. (1973a,b) and Kaiser (1973). The basic difference 
between the present work and that of Jones, et al . (1973a, b) is to use a 
propagator in place of U in equation (19) for D^(t). is the 
propagator in the partially averaged field, 

|3p(x,x') = f + §/ (x) • R(x-x') (93) 

A S^rther approximation is then made, namely, that in evaluating 
trajectories in £ , a guiding center approximation can be made (Kaiser 1973) 
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One consequence of this is that U (t) 0 U q (t) ^ = ^ +T w ^ete 0 is the 
phase angle. 

The propagator U(t), on the other hand, introduces the quantity 
Ap(t) through the relation 


U(T)p = P 0 (t) + Ap (t) 


(94) 


The approximation made by Jones et al. (1973a, b) that U p (T)0£j0+T 
is a fairly severe one. It is not difficult to show that such an 
approximation is equivalent to a Ap^ of the form (at jj, = 0) 


Ap x (T) = acos (0 + t) (95) 

Ap y (r) = asin (0 + t) (96) 

where we have chosen {3 to be in the z-direction ($ = e ) , and where 

z 

a is an arbitrary function of t and 0. However, Jones, et al. (1973a, b) 
make the additional approximation that 

(l-M,* S ) 1/a ~ a-H 0 S ) 1/2 = P oi (T) (97) 

where here 

U p (T)(Jl S (JL* 

u o (T)fx S M, o (98) 

But because the magnitude of the momentum is a constant, one 
must have 

(l-|j,* 2 ) = Cp 0 ^(t) + Ap . (T)] S (99) 

(we have dropped quantities of the form *<p ' (t) *p' (t)>) . Equations (97)- 

(99) together require that a = 0 (or equivalently, 6p (t) - 0). If 

one sets the term (t) = 0 in equation (88) , one finds the results 

-L 

plotted as open circles in figures (4) and (5). Now, in place of 
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equation (89) , one has 


€D (6p =0, |x=0) oc T] 3 (100) 

P- l 

We conclude then that the parameter dependence found by Jones, 
et al . (1973a, b) and Kaiser (1973, 1974) results from an inadequate 
approximation to the particles' momentum, p(t)» It remains to be shown 
whether more refined calculations of orbits in the partially averaged 
field would bring their theory into closer agreement with the work 
presented here. From the numerical simulations presently available, one 
cannot determine the parameter dependence of D with any precision. The 

r 

3 3 

published simulations are consistent with either an T| or behavior. 

Several important questions remain. First, the convergence 
properties of this theory have not been investigated. In the case of 
electrostatic turbulence, a convergence proof has been presented by 
Thompson and Benford (1973), but challenged by Orszag (1975). The 
electromagnetic case is considerably more complex. Perhaps comparison 
with numerical simulations is the most fruitful way to test this type 
of nonlinear theory. 

It is not difficult to generalize the results presented here for 
different correlation functions, including isotropic turbulence, and 
this work is in progress. The exponential correlation function is a 
convenient one to use for many reasons, but one should bear in mind that 
it has some well-known undesirable features. 

Also, the validity of the adiabatic approximation has not been 
investigated in any detail in the context of this resonance - broadening 
approach. It is known to fail at (j, = 0 in quasilinear theory, but is 
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expected to be a better approximation here because the time integral 
in equations (70) and (88) for D (<r) saturates quickly. These questions 
are presently being investigated. 
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VI CONCLUSIONS 


We have presented a theory for evaluating the pitch angle scattering 
coefficient for particles moving in magnetostatic turbulence. The theory 
is a perturbation expansion in the strength of the magnetic field fluctua- 
tions , which is equivalent to the weak coupling approximation of Kadomstev 
(1965). The pitch angle scattering coefficient was defined and evaluated 
assuming that the adiabatic approximation is valid. This theory is expected 
to be good for weak field fluctuations (Tf«l) and arbitrary particle 
rigidities. [However, for slab turbulence, the high energy regime, e«l, 
is not physically interesting.] We believe the derivation to be fairly 
rigorous with a minimum of ad hoc assumptions. 

The theory was evaluated in detail for the special case of a slab 
turbulence model with an exponential correlation function. For pr- 1 , we 
find that quasilinear theory provides a good approximation. At |jl=0, 
where quasilinear theory was expected to fail, this strong turbulence 
theory predicts significant scattering through 9CP — much of it due to 
mirroring forces. At p,= 0 , was computed in detail. All of the impor- 
tant nonlinear terms were retained. In that calculation (section III, 
above), in addition to the weak coupling and adiabatic approximations, 
it was assumed that terms proportional to P||(t) could be dropped. A 
test of the validity of that assumption will have to await more refined 
numerical simulations. In section IV, was computed for arbitrary p,. 

There terms in Ap (t) were omitted from the analysis; except the one 
term resulting from mirroring forces. The results of section IV 
evaluated at p *— 0 compare favorably to the more exact calculation of section 
III (cf. figs. 1-3) . 
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We found that eD^p-O) oc TJ 3 e . A direct comparison of these results 
with those of Kaiser (1973 5 1974) and Jones et al. ( 1973a ,b) showed that 
their approximate evaluation of their theory neglected the tern in Ae (t) 
that comes from mirroring. Our theory reduces to theirs (within factors 
of ~ 2) if one drops all terms in Ap x (t) in the analysis fcf. figs. 4 and 5). 
It remains to be shown whether a more accurate calculation of particle 
orbits in the partially averaged field of Jones, et al. (1973a ?b) and Kaiser 
(1973 j 1974) will bring their results into agreement with ours. The 
simulations presently available (Kaiser 1974) cannot distinguish between 
a pitch angle scattering coefficient proportional to T] 3 or eT| 3 . 

The success of this theory in matching the numerical simulations 
leads to the conclusion that the weak coupling approximation is a good 
one in this problem. The adiabatic approximation also appears valid, 
though more extensive checks of this are possible and are in progress. 


Our general conclusion is that can be accurately computed within 
weak coupling, strong turbulence theory, if a sufficiently good evalua*- 
tion of the particles’ momentum (and position) is included. We believe 
that this is the first calculation of to do this. Generalization 
to other correlation functions and other turbulence geometries 
isotropic) is in progress. 
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appendix a 


In this appendix we show that, for slab turbulence at p,=Oy- the 
.following terms are approximately zero 

<p! (v) -pi (\)>“ 0 (Al) 

<P^ (v) expik ||X|| (\)> ~ 0 (A2) 

In analogy to equation (42) we immediately have 

<Pi (v)*p!(\)> = -if JdrJds [[ U(v-t)p*Q' -C^t) J *^.(-s) *n’[x(X-s)] 
•p(X-s)> (A3) 

Using equation (47) for Q' (x), one has 

<£i(v)*£i(X)> = -T1 2 o J'dT o J’ds <p(v-T)*te , [x(v-T)]*n*C^(T+s).Q. 

P'[x (X-s)] (X-s)> (A4) 

If one writes the fluctuation fields as Fourier transforms, and performs 

the indicated matrix algebra, one finds 

* s v X » 3 

<£i(v)-£*(X)> = -p=r JdT J*ds JdkjjR(k|| ]cas(r+s) exp [-k^ (v-t,X-s) /2] 
V 2rr 0 0 -® 

<P|f (v-t )p,j (X-s)> (A5) 

where we have expanded the correlation and used the fact that at p,=0, 
< Pll( T ) >=:0 » In principle, it is possible to retain equation (A5), but 
to do. so would significantly complicate the analysis. However, this 
term is not expected to produce a significant contribution to Dj. (t) 
because the correlations in equation (A5) are in general evaluated at 
different times. In the context of the weak coupling approximation, such 
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nonlocal behavior is expected, to be small. Also, the correlation 
<p’ (v)Pjj (X)> is formally proportional to if , which implies that 
<pi’(v)«El'(X)> is ©(T] 4 ). For these reasons we arrive at equation (Al). 

We return to equation (A2), which can be written as 

<£i(v)expik||X|j (X)> = <2j.(v)expik||X|j (X)> (A6) 

Note that 

e' (v) s EU(v)-U(v)]p 1 (A7) 

From the identity (l6), we then have 

<£^(v)expikj|X|| (X) > = Uu(v-s)[X* u(s)-<£’U(s)>]p JJ 

expikj|X|| (\)> (A8) 

If one approximates U as U in equation (A8), then the second term is 
clearly of higher order in T) and one has 

<p^v)expik||X|| (X)> as -T) < q J’dr[p(v-T)£ , Uo (t)pjJ expik||X|j (X)> (A9) 

where we have used equation (36). If we rewrite expi[k||X|| (X)] as U(X)expi[k||X||] , 
and use equation (l6), equation (A9) becomes 

v _ X 

<p^(v)expikj|X|| (X)> as T] 2 < J\iT[ju(v-T).£ , bb ( T )Pxj] «fa s U(X-s)£’U(s) 
expik||X||> (Alo) 

Following Ben-Israel, et al^ (1975) we find, to lowest order in T|, 

v fi* _ a t i* -n ^ __ 

<p^’(v)expik||X|| (X)> a* -iT| s J*dT < [j_ U(v-t)j>Q'6 l (t)Jj Jds U(X-s) 

(All) 
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V /V f I 

JdT J\is R(k||)<expik||[x||(v-T)-X||(X-s)]sk||P||(v-T)£ A (x- s ) 

)> (A12) 

In the analysis presented here, we will drop any terms in £]|(t). 

This is equivalent to ignoring nonlocal correlations of the form <£* (t)*' (\)>. 
we expect that in the context of the weak coupling approxima- 
tion, nonlocal correlations between momentum and position fluctuations 
should be less important than autocorrelations of fluctuation in posi- 
tion evaluated at equal times [ e.g. Pi ( t,t)= <x f (t)x’ (t )>] . For a more 
comprehensive discussion of these various correlations, and their role 
in determining D^(t), the reader is referred to Ben-Israel, et al, (1975). 
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APPENDIX B 


In this appendix we derive equation (55) for Ap, (t) . The result 
of taking the Fourier transform of the fluctuation fields in equation 
(5*0 is that 


2 T T 

Ap (t)^ — 3— — PdX Pdv fd 3 k J , d 3 k , <exp[ik«x(T-v)]p(T-v) *0’ (k!) ‘^(t-v) • 
(2tt) 3 °. X 

p-- + n» 5 *(T)»vJ p(T-X)»fl , (k , )* 6 *(\)exp[ik , «x(T-X)] > (Bl) 

dp “ J “ 1 


where we have used the approximation 


U(t) —£• — U(t)sC(t)-[^j- + n.^ t (T)-v]u(T) (B2) 

ap apb(t) dp 

If one substitutes the identity (47) for O’ (k) and uses equations (22) 
and (44), it follows, after some straightforward algebra, that 


2 T T « . , - , • 

i B±( T ) = -^= — JdX x J’dvJ’dkj|R(k|| )<exp[ iky xy (v-X)] jp|| (t-v)C(v-X-t) 


TV Pll( T-X ) + ik||P|| (t-X)V- rj| (t-X)]-p x (t-v)*C^ (X) 

dpi dp x 

TV PH (t-\) + ik||P|I (T-X)( V rj|(T-X)-(T-v))| > (B3) 


Again, we drop terms containing pjj (t ) . Equation (55) follows immediately 
when one uses the approximation (A 2). 
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APPENDIX n 

Prom equation (4o) for A£(t) we have 

T T 

Apu(t) = T| a JdX J\iv <U(t-v)p»Q* - r U(v-X)p*f?»|^> 
o X ^p 

Following the derivation of Axj|(T)given in §IV, one has 

2 J Y ® 

AP|| (T)=-^zr“ n Jd X < fdvJ'dk||R(k|| )exp{ik|j[|J.(v-X)]-Ax|| (v-X)] 
\/2tt 0 A 

£HT'(v-X )[2p,-ik|| (l-|j, a ) ( v-X )] 
which becomes 

Apjj (t) = ^rr JdX J’dk||R(k|| )exp[-k|D 1 (X )] (t-X) 

0 T 

{2\i iar[k|| (|iX-Ax|| (X))] + (l-y® )kjjX sin [ k|j (pA-Ax|| (X))]} 
Equation (56) follows immediately. 


( 01 ) 

*iDi(v-X)/2} 

( 02 ) 

( 03 ) 
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APPENDIX D 


In the discussion below we outline the transformation of equations 
(70)-(72) and (63)-(64) to a system of integro-differential equations 
that can be solved numerically. Equation (70) can be easily converted 
into a differential equation 

Dp i ( T )=^ ]D ( j L ( T )- - ^ : ^rg( T )[ har ( T ) + Ap x (t)] (Dl) 

Differentiating equation (71) yields 

T T T 

d i(t)=-§ 2 — J<3XwXg(X)(T 2 -XT)+-£^ JdX Jdvg(v-X) (2 t-v-\) 

V2tt 0 \J2tt ° o 

CA£ x (^) , P J£) (X) + ^A£ x (v)-A£ 1 (X)] (D2) 

Now define 

fi (t) = g(T) cost (D3) 

f a (t) = tg(T) cos t (d4) 

£ 3 (t) = JdX J'dvg(v-X)AE x ( v )-[pj 0 (X)*tiA£ x (X)] (D5) 

0 o ' 

f 4 (T) = JdX J 1 dvg(v-X)(v+X)A£ JL (v)-[p. (X)+iA£ 1 (X)] (P6) 

0 0 ^ 

with f^O) =0, i = 1-4. Therefore 

Dj. (r) = -^£- [T®fi (T)-Tf a (T)+2Tf 3 (T)-f 4 (T)] ( D 7) 

\[2TT 

Now cons ider 

T A t \ 

f 3 (t) = JdXg(X)A£ ,(X).p Jk) (T-X)+ Jdvg(v-T)A Ex (v)*p. (t) 
o o 

T 

+ JdXg(T-X)A£ x (X)*£ x (T) ( d8 ) 

o 
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and define 


ffi(T) = g(T)sin T (d.9) 

fe ( T ) = Tg(T)slnT, |_f B (0) = f 6 (0) = 0] (DIO) 

T 

<3i (t) s J v dvg(T-v)Ap lX (v) (Dll) 

T 

1s(T)= o J‘dvg(T-v)Ap iy (v) (D12) 

T 

<la(T)= Xdvvg(T-v)Ap >U£ (v) (D13) 

T 

q4( T ) s o J’ d v v g( T -v)AP iLy (v) (Dl4) 

so that 

• - 

f 3 (T) = fi (t)[Ap j:jc (t)(^5't + Ap X y(r)sinT ] + 
f s (T)[Ap xx (T)sin t - APj^tVkjt] + 

1 1 (t)[«U't + Ap lx (r)]ti3(T)[sin t + Ap^t)] (D15) 

Similarly, 

f 4 (T) = [2Tfi (T)-f 2 (T)][Ap xx (T)%5 > T + Ap^T )sin t] + 
[2Tf 5 (T)-f s (T)][Ap ix (T)sln T - Ap jy (T>WT] + 
[Tqi(T)-tq 3 (T)][fl«*T +' ; APix(t)] + [Tq 2 (t)*U (t) ][ sln T+Ap_^.(t)] 

(Dl6) 

In a completely analogous manner one has 
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4 Pi x (T) = -pr- tafTCfxCOU + t- AP|| (t) ) + Tf s (T)-f 6 (T)+f 7 .(T) + 

\[2tt om< H -=0 

+ f3o(T)+f»(T)] 

-Sin T[-f s (T)(l + Apn (t) )^ =0 + Tf l(T)“f , 3 (T)+f 8 (T) + 

+fg(T)+f 11 (T)] + q 1 (T)[l + ^ AP||(t)]^ =0 3 (D17) 

2 

Ap jy (T)=-^r {^T[f 5 (T)(l + Ap||(T)) Mi=0 -Tf X (T)+f s (T)-f 8 (T) 

-fg(lO)-f n (T)] 

-sin t{* f*! ( t ) ( 1 + Ap b (t)) + Tfg (t ) ~ f 6 (t) +f *7 (t ) 

+f 10 (T)+f te (t)] + qg (t)[1 + AP||(T)]} (D18) 

oM- 

where 


^7 (T) = f* B ( T )|^ AP||(t) 

H=0 

(D 19 ) 

h(T) = f l ( T )|^ Apjj (t) 

(JL=0 

(D20) 

fg Ct) = qj, (T^T-qg (r)sin t 

(D21) 

F]o(t) = qg (T^T+q! (T)sin t 

(D22 ) 

%(t) = fg(T)[|ji A P || (t )]^_ 0 

(D23) 

f 3s( T ) = A P ||(T)r =0 

(D24) 


Finally 
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where 


Aph(t)|. - [Tf 13 (T)-f u (t)] 


i=0 \[2tt 


( 52 5 ) 


foa(T) = iMTT Cs(t) + Th(T)[T - AX|| (t) 1^ =q } 
*1a(t) = Tf^T) 


(D26) 

(D27) 


and 


t?'^U = -^ 45 » (T) l=c 


(D28) 


The complete system of equations (Dl) , (D3)? (D4), (D7) ? (D9) a-* 1 * 3 - 
(D28) can he numerically integrated using a standard Adams -Bashf or th 
predictor - corrector method. 


44 



APPENDIX E 


The system of equations derived in §IV can be converted to the 
following set of differential equations. The derivation is straight 


forward. 

= ^=r (1-M' 2 )G(t)[1+6 Px (t)] 5 D (0)=0 (El) 

^ V2n P 

Di (t) = T[Te 1 (r)-e 3 (T )] 3 D : (0)=0 (E2) 

V2tt 

©i (t) = G(t) (E3) 

©s (t) = tG(t) (e4) 

Ax|| (t)=-^- C2M,[Te 1 ( T )-e 3 (T)] + (l-|i 3 ) [Te 3 (t) -e 4 (t) ] } (E5) 

\[2tt 

k (t) = H(t) (E6) 

©4 (t) = tH(t ) (E7) 


G(t) and H(t)are given by equations (78) and (80). 6p x (t) follows from 
equation (86). 

% 
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FIGURE CAPTIONS 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


The dimensionless pitch angle diffusion coefficient 6D 

P 

is plotted as a function of |j,. The solid square represents 
the solution of the set of equations derived in §111, While 
the solid line follows from the system of equations 
discussed in §IV. The parameter values are 6=1, /2T|=0.1. 

6D vs. )jl for 6=2, /2T)=0.15. The notation follows that 

r 

of figure 1. 


€D vs. for 6=1, /2'T1=G.3. The notation follows that of 

P “ 

figure 1. 

€D vs. 6 at u=0 for /2T|=0.3. When the momentum deviation 
P 

terms are retained (Ap(t)^O), €D ((j,=0) is approximately 

— P 

proportional to 6 (solid circles)* However, when these 


momentum terms are dropped in equation 88) one 

finds that €D (p<=0) is approximately independent of € 
P 

(open circles), thus recovering the results of Jones, 


et al . (1973a, b). 

€D (|i=0) vs. T| for €=5. €D (€=0) is approximately 
P P 

3 

proportional to T] whether the momentum deviation terms 
are kept (solid circles) of dropped (open circles). 


48 



JE_} PRESENT WORK 

JONES, ET AL. ( 1973 a,b ) 

QUASI LINEAR THEORY 


<g) VOLK (1973) 

X OWENS (1974) 

NUMERICAL SIMULATIONS (KAISER 1974) 





5 X I0“ 3 


4X10* 


€ UfJL 


3X 10* 


2XI0“ 3 


IX 10* 








€ =1 



€ Dfi ifi - 0 ) 




D/I ( LL =0 ) 



OF POOR QUALITY 



